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ABSTRACT 

Context. Empirical evidence suggests a tantalising but unproven link between various indicators of solar activity and the barycentric 
motion of the Sun. The latter is exemplified by transitions between regular and more disordered motion modulated by the motions of 
the giant planets, and rare periods of retrograde motion with negative orbital angular momentum. An examination of the barycentric 
motion of exoplanet host stars, and their stellar activity cycles, has the potential of proving or disproving the Sun's motion as an 
underlying factor in the complex patterns of short- and long-term solar variability indices, by establishing whether such correlations 
exist in other planetary systems. In either case, these studies may lead to further insight into the nature of the solar dynamo. 
Aims. Some 40 multiple exoplanet systems are now known, all with reasonably accurate orbital elements. The forms and dynamical 
functions of the barycentric motion of their host stars are examined. These results can be compared with long-term activity indicators 
of exoplanet host stars, as they become available, to examine whether the correlations claimed for the Sun also exist in other systems. 
Methods. Published orbital elements of multiple exoplanetary systems are used to examine their host star barycentric motions. For 
each system, we determine analytically the orbital angular momentum of the host star, and its rate of change. 

Results. A variety of complex patterns of barycentric motions of exoplanet host stars is demonstrated, depending on the number, 
masses and orbits of the planets. Each of the behavioural types proposed to correlate with solar activity are also evident in exoplanet 
host stars: repetitive patterns influenced by massive multiple planets, epochs of rapid change in orbital angular momentum, and 
intervals of negative orbital angular momentum. 

Conclusions. The study provides the basis for independent investigations of the widely-studied but unproven suggestion that the Sun's 
motion is somehow linked to various indicators of solar activity. We show that, because of the nature of their barycentric motions, the 
host stars HD 168443 and HD 74156 offer particularly powerful tests of this hypothesis. 
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1. Introduction 

The centre of mass of the solar system, the barycentre, moves 
through space controlled by the Galaxy's gravitational potential. 
For solar system dynamical studies, the barycentre is the origin 
of a local quasi-inertial frame. The Sun, the planets, and all other 
gravitating matter orbit this centre of mass, resulting in a contin- 
uous displacement of the Sun's position relative to the barycen- 
tre. This displacement is determinable, with adequate precision 
for the objectives of this study, by Newtonian dynamics. 

As noted in the first discussions of the Sun's orbital motion 



about the barycentre by Newton (as quoted by Cajori 1934 1, the 
motion of the Sun is rather complex 'since that centre of gravity 
is continually at rest, the Sun, according to the various positions 
of the planets, must continuously move every way, but will never 
recede far from that centre.' 

Jose| ( 1965) used the improved numerical integration of the 



outer solar system planets over about 400 years by |Eckert et al.| 



( 1951| l to determine the path of the Sun with respect to the solar 
system barycentre more precisely than had been possible pre- 
viously. He claimed a 178.7 yr periodicity in its principal fea- 
tures, and suggested correlations of this, and the rate of change 
of the Sun's orbital angular momentum, with both the 1 1-year 
sunspot cycle, and the 22-year magnetic dipole inversion ( ]Hale| 



1924) and sunspot polarity cycles. Juckett (2000) criticised the 
latter claim on the basis of Jose's arbitrary phase adjustments. 
Whether these specific correlations exist or not, by drawing at- 
tention to various periodicities in the Sun's barycentric motion 
|Jose| ( |1965] l nevertheless stimulated further research into its pos- 
sible link with solar activity. 

A substantial body of work has since appeared on the pur- 
ported link between the Sun's motion and aspects of its long- 
term activity. With only one exemplar, the subject has reached 
something of an impasse. We show that exoplanet systems pro- 
vide numerous other examples to allow further progress. 

2. The case of the solar system 

In the simplest case of a single orbiting planet, the star and planet 
displacements from the system barycentre are related by 



r*(0 



r P (0 ■ 



(1) 



where M p and M+ are the planet and star masses in common 
units, and r p and r* are the positions of the planet and star with 
respect to the barycentre, again in common units. Both planet 
and star move in elliptical orbits in the same plane, about the 
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common centre of mass, with semi-major axes which differ by 
the ratio M p /M*, and with orientations (specified by their argu- 
ment of pericentre, to) which differ by n rad. 

In a multiple planetary system in which the planets are grav- 
itationally non-interacting, and in which the individual and total 
planet masses are small with respect to the mass of the star, the 
motion of the star about the system barycentre can be approxi- 
mated by the linear superposition of the reflex motions due to the 
Keplerian orbit of each individual planet around that star-planet 
barycentre. If the planets have periods or close approaches such 
that they are dynamically interacting, the orbits must then be de- 
termined, at least over the longer term, by N-body integration. 
This is a complication that we ignore in this study. 

Equation 1 shows that low-mass, short-period planets natu- 
rally have a smaller effect on the Sun's (or star's) displacement 
compared with more massive planets, and particularly those with 
large semi-major axes (i.e. long orbital periods). 

In the case of the solar system, the Sun's motion under the 
influence of the four giant planets (Jupiter, Saturn, Uranus, and 
Neptune) is illustrated in Figure [2^ for the time interval 1970- 
2030. Equivalent plots, with which our own determinations are 
consistent, can be found in , e.g.,|Jose|( 1965 Figure 1) over the 
interval 1834-2013, and in |Juckett| ( |2000| Figure 2) over the in- 
tervals 1 8 1 5-65 and 1 9 1 0-60. 

The Sun's barycentric motion is evidently most strongly in- 
fluenced by the motions of these four giant planets. Patterns 
in rotationally invariant quantities, such as successive close ap- 
proaches between the Sun and the system barycentre, termed 



We expand briefly on two of these themes to set the context. 
One of the specific curiosities of the Sun's motion noted by |Jose| 
(1965) is evident in Figure [2^. Around 1990, and before that in 
1811 and 1632, the Sun had a retrograde motion relative to the 
barycentre, i.e. its angular momentum with respect to the centre 
of mass was negative ( |Jose|[T9 65). The next such retrograde Sun 
event will occur around 2169. In a multi-body system in which 
total angular momentum is conserved, such behaviour of one 
particular body is, again, not in itself surprising. The result is of 
more particular interest because Javaraiah ( 2005 ) has argued that 
epochs violating the Gnevishev-Ohl 'sunspot rule', which states 
that the sum of sunspot numbers over an odd-numbered cycle 



'peribacs' by |Fairbridge & Shirley (T987J, recur, for example, 
at mean intervals of 19.86 yr, corresponding to the synodic peri- 
odic of Jupiter and Saturn. 



2.1. Relation to solar activity 

Solar axial rotation plays a fundamental role in the two main hy- 
potheses for a mechanism underlying the solar cycle: attributed 
to a turbulent dynamo operating in or below the convection en- 
velope, or to a large-scale oscillation superimposed on a fossil 
magnetic field in the radiative core. The precise nature of the 
dynamo, and many of the details of the associated solar activ- 
ity (such as the details of the sunspot cycles, or the reasons for 
the prolonged Maunder-type solar minima) remain unexplained 
(e.g. Charbonneau, 2005 1 Choudhuri, 2007; Weiss & Thompson 
2009 ; |Jones et al.||2010[ ), although certain features may arise nat- 
urally in some models (e.g. |Charbonneau et al. 2007[ ). 

Empirical investigations have long pointed to a link between 
the Sun's barycentric motion and various solar variability indices 
(e.g. |WoIf| [TB591 |Brown[ [T900] |Schuster| [T9TT] |FerrTsl [T969l >. 

For the reasons noted below, such a connection has not at- 
tracted much serious attention in theories developed to under- 
stand the solar dynamo. Nevertheless, acceleration in the Sun's 
motion, or the change of its orbital angular momentum, have 
been linked empirically to phenomena such as the Wolf sunspot 
number counts (|Wood & Wood |1965| l, climatic changes (Morth 



exceeds that of its preceding even-numbered cycle ( Gnevishev & 
OhTj [T9481 IKomitov & Boney] |200T] [Charbonneau et al.| |2007t 
Nagovitsyn et al. , 2009 ), are close to these intervals of the Sun's 
retrograde orbital motion. 

|Charvatov a (2000) drew attention to periods of the Sun's or- 
bit which are characterised by more ordered 'trefoil'-type mo- 
tion dominated by Jupiter and Saturn, and which recur with the 
periodicity of 178.7 years. More 'disordered' motion appears 
between these periods (their Figure 2). |Charvatova] ( |2000| ) as- 
sociated these intervals with the prolonged minima of the solar 
activity cycle, viz. the Wolf (1280-1350), Sporer (1430-1520), 
Maunder (1645-1715), and Dalton (1790-1820) minima. 

2.2. The question of a causal link 

The implication of any relation between solar motion and solar 
activity, should it exist, is that there is some mechanism related 
to its motion which, in some way, alters or modulates material 
flows and/or magnetic fields within the rotating and revolving 
Sun. In Newtonian gravity, only planetary tidal forces (oc r~ 3 ) act 
externally to stress or distort the figure of the Sun, and these have 
been excluded as the underlying causal factor because of their 
small amplitudes ( |Okal & Anderson] 1975 |Smythe &~E ddy 



& Schlamminger| 1979 1, the 80-90-yr secular Gleissberg cy 
cles ( |Landscheidt| |1981| |1999|l, the prolonged Maunder-type 



solar minima (Fairbridge & Shirleyj |1987| |Charvatova| 1990 



2000 ), short-term variations in solar luminosity ( |Sperber et al 
1990 1, sunspot extrema (Landscheidt, 1999), the 2400-yr cy- 



2000), hemi- 



cle seen in 14 C tree-ring proxies ^Charvatova 
spheric sunspot asymmetry ( Juckett 20 00[), tor sional oscillations 
in long-term sunspot clustering (Juckett, 2003 ), and violations of 
the Gnevishev-Ohl sunspot rule (Javaraiah 2005 [ ). 



1977). In metric theories, that the Sun is in gravitational free 
fall has also been used to argue that there should be no causal 
link between its barycentric motion and any type of solar activ- 
ity indicator ( |Shirley||2006] ). 

Attempts to identify a spin-orbit coupling mechanism be- 
tween the solar axial rotation and orbital revolution neverthe- 
less continue. Thus Zaqarashvili (1997]) has suggested that the 
Sun's predominantly elliptical motion creates a periodic velocity 
shear on the weak poloidal magnetic field in its interior. |Juckett| 
(2000 ) has instead proposed that the orientation of the Sun's ro- 
tation axis (inclined by ~7.°25 with respect to the ecliptic, and by 
~6.°25 with respect to the invariant plane, viz. the plane defined 
by the total angular momentum of the planets and the Sun) facil- 
itates non-linear frequency mixing of planetary-induced forcing 
oscillations, and induces differential flows to conserve angular 
momentum. 

Although elements of each explanation have been criticised 
by |Shirley ( |2006[ ), the latter author nevertheless argues that a 
clear physical picture, properly linking the Sun's rotation and its 
revolution about the barycentre, remains elusive. 

We can summarise by stating that an empirical connection 
between solar activity and the Sun's motion appears intriguing, 
but remains unproven and unexplained. Perhaps like the Titius- 
Bode law, which has been largely attributed to phenomenolog- 
ical patterns which can be replicated by broad families of two- 
parameter geometrical progressions (Hayes & Tremaine 1998; 
Murray & Dermott 2000, Section 1.5), it is purely specious. But 
it has been extensively hypothesised, and an independent oppor- 
tunity for testing it would be of interest. 
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Fig. 1. Multiple (Doppler) planet systems, ordered by host star 
mass (indicated at left with size proportional to M*, ranging 
from 0.31M G for GJ 581 to 2.8M Q for BD +20 2457). Each 
planet in the system is shown to the right, with sizes propor- 
tional to logMp (ranging from about 0.01 - 20Mj). Horizontal 
bars through the planets indicate maximum and minimum star- 
planet distance based on their eccentricities. 



3. The barycentric motion of exoplanet host stars 

Multiple exoplanet systems, those with two or more known or- 
biting planets, provide an independent possibility of testing the 
hypothesis that the Sun's barycentric motion is somehow linked 
to solar activity. To do this, specific features of their host star 
barycentric motions (as determined through Doppler measure- 
ments in the radial direction) can be compared with their own 
stellar activity indicators. This study is concerned with the for- 
mer aspect. 

Of nearly 500 known exoplanets, 41 multiple systems dis- 
covered by radial velocity measurements are currently known. 
Of these, 30 systems have two known planets, while 1 1 systems 
have three or more. This sample of Doppler-detected multiple 



Table 1. Barycentric host star motions, ordered by increasing 
| dL,/df | max . Maximum and minimum values of the star's orbital 
angular momentum are given for the time interval 2000-2100. 
lM AU 2 d-' = 5.15 x lO^kgmV; lM AU 2 d- 2 = 5.96 x 
10 42 kgm 2 s- 2 . 
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systems is complemented by the two multiple systems discov- 
ered from photometric transits, CoRoT-7 and HAT-P-13. 

To determine the host star barycentric motions, we have 
used the orbital elements for the multiple systems compiled 
at www.exoplanets.org. The characteristics of the Doppler- 
detected systems according to stellar mass, planetary mass, or- 
bital separation and eccentricity, are shown schematically (but to 
scale) in Figure[T] The eccentricity, e, semi-major axis, a, and ar- 
gument of pericentre, to, define the individual orbit geometries. 
The time of periastron passage, f p , and orbital period, P, addi- 
tionally determine the position of the planet along its orbit. 

We assume lower limits for the planetary masses, i.e. M p sin i 
with sin i = 1 (giving lower limits to the stellar displacements), 
and further that all planets in a given system are coplanar. We 
ignore dynamical planetary interactions known to be significant 
in some cases, and quadrupole or relativistic precession, as irrel- 
evant for the qualitative orbital properties over relatively short 
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time intervals. For the Sun, only the contributions of Jupiter, 
Saturn, Uranus and Neptune are considered. Their eccentricities 
are included, although they have little additional effect. 



four orders of magnitude larger than the negative excursion of 
the Sun (not present in the time interval of Table [TJ. Stars with 
the next most negative extrema are HD 74156 and HD 217107. 



3.1. Orbital angular momentum 



4.1. Orbital motions 



In the orbit plane, the contribution of the star's displacement with 
respect to the system barycentre r, its velocity v, and its accel- 
eration v, are calculated for each planet i, under the assumption 
Mi «: M+, as 



Mi di 
M tn , 



cos to(e - cos E) + Vl - e 2 sin to sin E 

sin to(e — cos E) — Vl - e 2 cos co sin E 




(2) 



2nMj a; 



Vj = - 



M tot P(l -ecosF) 
An 1 Mi a -, 



cos to sin E + Vl - e 2 sinwcosii' 
sin oj sin E 

_ VTT 

<? z COSWCOS/i (3) 





M tot P 2 (l -ecos£) 3 



cos to(e - cos E) + Vl - e 2 sin u> sin E 

sin co(e - cos E) — Vl - e 2 cos w sin E 




(4) 



where M, are the individual planet masses, M tot = M* + £ M,, 
and a, are the semi-major axes of the planet's orbit around 
the star (not the barycentre). In these equations, the eccentric 
anomaly E at time t is derived iteratively from Kepler's equa- 
tion, 



E - e sin E 



2n 



<t ~ f P ) ■ 



(5) 



The total angular momentum L of the star with respect to the 
system barycentre, and its time derivative (identically zero for a 
single planet system), are then given by 



L 

dL 



= M+ r x v , 
= M* r x v , 



(6) 
(7) 



where r, v, and v are the sums over the individual contributions, 
from which we derive the components perpendicular to the com- 
mon orbital plane 



L z = M+(r x v y - r y v x ) , 
dL z 

— = MJr x v y - r y v x ) . 
at 



(8) 
(9) 



In the accompanying plots, the xy-projection corresponds to 
the orbital plane, and the x-axis points from the barycentre to- 
wards the common ascending node. In the case of the solar sys- 
tem, the xy-plane essentially corresponds to the ecliptic, and the 
x-axis points from the barycentre towards the equinox. 

4. Results 

Table[T]summarises the results from the dynamical analyses, or- 
dered by increasing | dL z /df | max determined over 100 years from 
2000-2100. All but two of the stars have | dL./df | max exceeding 
that of the Sun. HD 168443 has the largest dL,/df, more than 
five orders of magnitude larger than that of the Sun. 

The majority of stars, 26 out of 4 1 , also experience periods of 
negative orbital angular momentum, L z . HD 168443 is further- 
more the system with the smallest (most negative) value, some 



Figures |2]-[3] show a selection of the resulting host star barycen- 
tric paths, including that of the Sun. The behaviour can be com- 
pared with the architecture of the systems shown in Figure [T] 

The path of the Sun (Figure|2^) replicates that of earlier work 
(e.g. |Jose| |1965| |Juckett 2000| l, with maximum displacements 
exceeding 1R Q . Our determinations of L z and dL z /dt also closely 
match those of |Jose| ( [l965 their Figure 2b) and Javaraiah (2005 



their Figure 1). An excursion to negative L- occurs around 1990. 

A number of systems, including those dominated by mul- 
tiple planets with a < 0.2 - 0.3 AU, result in only small per- 
turbations from a predominantly elliptical stellar path, even for 
some triple or quadruple planetary systems. Examples (not il- 
lustrated) are GJ 876, BD -08 2823, HD 11964, HD 38529, 
HD 47186, HD 128311, HD 187123, HD 190360, HD 215497 
and HD 217107. The five-planet system 55 Cnc is qualita- 
tively similar, displaying essentially circular motion dominated 
by the most massive outer planet (4Mj at 6 AU), but with low- 
amplitude high-frequency perturbations due to the second most 
massive planet (0.8Mj at 0.1 AU). 

Slightly more complex motion is seen in 47 UMa, HD 9446, 
HD 45364, HD 60532, HD 73526, HD 82943, HD 202206, as 
well as in the transiting multiple CoRoT-7. 

BD +20 2457 (Figure |2j>) is the most massive star in the 
present sample at ~2.8M . With two massive planets (23Mj at 
1 .5 AU and 1 3Mj at 2 AU), it displays a barycentric motion qual- 
itatively most similar to that of the Sun. 

HD 168443 (Figure|2};) has two massive planets at 0.3-3 AU 
(8 - 18Mj), with a high eccentricity e ~ 0.5 of the inner planet. 
The star shows the most extreme values of negative L z and 
largest dL z /df, despite its relatively simple orbit. These excur- 
sions are regular, with a periodicity of ~58 d given by the period 
of the innermost planet. Similar motion is seen in HD 74156, 
HD 147018, v And, and in the wide-transiting HAT-P-13. 

For fi Ara (Figure |2jl), a periodic looping motion is domi- 
nated by the two most massive outer planets (1.9Mj at 5.3 AU 
and 1.7Mj at 1.5 AU). The third most massive planet (0.5Mj 
at 0.9 AU) imposes more abrupt changes in orbital angu- 
lar momentum on the underlying two-planet induced motion. 
Similar effects are seen in HD 12661, HD 108874, HD 134987, 
HD 169830, HD 181433 and HD 183263. In HD 181433, the ec- 
centricities of the two most massive outer planets result in a host 
star motion systematically displaced from the system barycentre. 

For 61 Vir (Figure |2£), the three planets are of low mass 
(0.01-0.07 Mj), and lie at small separation (0.05-0.5 AU). While 
the motion is complex, the barycentric displacements and angu- 
lar momentum changes are all small. GJ 581 is similar. 

HD 37124 (Figure|2p is a triple system with almost identical 
planetary masses in the range 0.6-0.7Mj. A simple looping mo- 
tion like that of fi Ara results from the two most massive planets 
at 0.5 and 3.2 AU. But this is transformed into a complex spi- 
ralling pattern in the presence of the third planet at 1 .7 AU, aris- 
ing from the various planetary orbital periods of 154, 844, and 
2300 d. Somewhat similarly complex but low-amplitude struc- 
ture is seen in a number of other systems, including HIP 14810, 
HD 40307, HD 69830, and HD 155358 (Figure [5). 

On a historical note, a curve similar to the symmetric 
part of Figure |2jl, illustrating the orbit of Mars viewed from 
Earth, appeared in Kepler's 1609 Astronomia Nova. These 
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Fig. 2. Barycentric motion of the host star for a selection of representative multiple exoplanet systems. Main plots (central two 
columns) show the orbit over the indicated time interval in a reference frame with the system barycentre at the origin, with abscissae 
and ordinates in AU (the solar diameter is R Q = 6.96 x 10 s m, or 0.00465 AU). To the outer side of each orbital sequence, plots 
show the orbital angular momentum, L- (upper), and dL./df (lower) for the same time interval, in the units as given in Table[T] For 
the systems shown, star masses lie in the range 0.85 - 1.15M Q except for BD +202457 which is ~2.8M Q . Orbital parameters were 



taken as follows: the Sun: JSeidelm ann (2005); BD +202457: Niedzielski eta!] pD09) ; HD 168443: |Wnght et al.| ( |2009) ; fi Ara: 
Pepe et al.| ( |2007) ; 61 Vir:|Vogt et al.|(|2010); and HD 37124:|\fegtet al.|(|2005). 



families of curves have been named 'planet mandalas' by 
N. Michelsen & M. Pottinger, after the Sanskrit for circle, 
and appear as such in the Wolfram Demonstrations Project, 
http://demonstrations.wolfram.com/SolarSystemMandalas/ 

5. Correlation with stellar activity 

A detailed investigation of any correlation between these mo- 
tions and stellar activity lies beyond the scope of this paper. But 
we note that observable diagnostics of stellar activity are avail- 
able in principle, and are becoming so in practice. Amongst these 



are radial velocity 'jitter' -type measurements probing stellar at- 
mospheric inhomogeneities such as spots and plages modulated 
by stellar rotation (e.g.|Saar & Donahue|ll997|ISaar et al.[ |1998[ 
|Jenkins et al.] [2006] |Desort et al.| |2007l |Saar| |2009) ; chromo 



spheric activity measured through the Call H and K lines (e.g. 
|Saar & Cuntz}|2001| [Wright et al.||2004), sometimes attributable 
to tidal or magnetic interactions (|Cuntz et al.| |2000[l, with some 



monitoring programmes extending back to 1966 ( [Baliunas et al. 
1995|); star spots observed photometrica lly during transits (e.g 



Silva-Valio 2008 , Dittmann et al. 2009); and changes in photo 



spheric radius due to varying magnetic activity (Loeb 2009 1. 
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Fig.3. As FigureU showing further examples of barycentric orbits for HIP 14810, HD 40307, HD 69830, and HD 155358. 
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2005 



p - 0.048 AU (| Butler et aO [2006) . The star has a 
~Tof about 0.5Mj ( |Shkolnik 
3xl0" 4 T 



a et a 



2008 ). It has a weak magnetic field of 1 
| |2007| l, which underwent a global magnetic polar 



ity reversal between 2006-07 (|Donati et al. 2008[), and again 



between 2007-08 (Fares et al. 



dominantly poloidal and toroida 



2009 ), switching between pre- 



forms. Observations suggest a 



magnetic cycle of ~2 yr, much shorter than the Sun's 22 yr. 

[Fares et aL] ( |2009] > hypothesise that this is governed by the 
close-orbiting planet, through tidal synchronisation of the outer 
convective envelope enhancing the shear at the tachocline. 

From our study HD 168443, with the innermost planet at 
a p ~ 0.3 AU, offers particularly interesting prospects for investi- 
gation as to whether much more distant planets can affect stellar 
activity through the star's reflex motion. Present physical under- 
standing says that they cannot. But its dL z /df, with a periodicity 
of 58 d, exceeds by more than five orders of magnitude that of 
the Sun. If orbital angular momentum variations play a role, its 
effects should be visible in this system. 

HD 74156 displays the next most negative orbital angu- 
lar momentum intervals. Qualitatively similar to HD 168443 
(Figure [2};) these negative excursions extend over 27 or so 51-d 
periods (3.8 yr), and are then absent for the subsequent 20 or so 
periods (2.8 yr). Its large mass, 1.24M , further implies a shal- 
low convective envelope, perhaps with greater susceptibility to 
shear effects. If retrograde angular momentum plays a role, these 
specific modulations will offer a defining diagnostic. 

Such periodicities observed in their stellar activity would 
support the unorthodox picture that the barycentric motion of 
the Sun somehow affects long-term activity cycles. Otherwise, it 
would imply that any link between the two is spurious. 

Much larger orbital angular momentum is experienced for 
stars in wide binary orbits, including exoplanets in circumpri- 
mary (S-type) orbits such as that around 16 Cyg B, where the 
stellar binary has P or b ~ 18 000 yr. Heuristically, it is difficult to 
imagine that L z alone can have any relevance. 

6. Conclusions 

Our study demonstrates that a variety of complex barycentric 
motions exists for exoplanet host stars. Behaviour cited as being 
correlated with the Sun's activity, for example intervals of more 
disordered motion, large changes in orbital angular momentum, 
and intervals of negative orbital angular momentum, are com- 
mon - but more extreme - in exoplanet systems. 

Accompanied by detailed studies of the associated stellar ac- 
tivity, these systems offer an independent opportunity to corrob- 



orate the hypothesised link between the Sun's barycentric mo- 
tion and the many manifestations of solar activity. In particular, 
they offer the possibility of independently testing any theories of 
spin-orbit coupling which are advanced in the case of the Sun. 
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